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We have extended the validity of the correlation between the surface 3d 
core-level shift (SCLS) and the surface d band shift (SDBS) to the entire Ad 
transition metal series and to the neighboring elements Sr and Ag via accurate 
first-principles calculations. We find that the correlation is quasilinear and 
robust with respect to the differencies both between initial and final-state 
calculations of the SCLS's and two distinct measures of the SDBS's. We 
show that despite the complex spatial dependence of the surface potential 
shift (SPS) and the location of the 3d and Ad orbitals in different regions 
of space, the correlation exists because the sampling of the SPS by the 3d 
and Ad orbitals remains similar. We show further that the sign change of 
the SCLS's across the transition series does indeed arise from the d band- 
narrowing mechanism previously proposed. However, while in the heavier 
transition metals the predicted increase of d electrons in the surface layer 
relative to the bulk arises primarily from transfers from s and p states to d 
states within the surface layer, in the lighter transition metals the predicted 
decrease of surface d electrons arises primarily from flow out into the vacuum. 
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I. INTRODUCTION 



The surface electronic structure of a metal plays a central role in the surface chemistry 
that takes place upon it. That is, the surface electronic structure of a metal is central 
to its chemical reactivity. It has long been understood that in the transition metals in 
particular, the d electrons dominate their chemistry. More recently, computational studies 
of chemisorption and dissociation on various metals and alloys have clarified how they do 
so (T|-§[. Thus, one can state with confidence that the surface chemical reactivity of a 
transition metal depends strongly on the response of its surface d electrons to the external 
perturbations imposed by an atom or molecule with which it interacts. A generic measure 
of its reactivity should then be provided by their response to a generic perturbation. 

One very simple generic perturbation is provided by the difference in environment be- 
tween surface and bulk. Thus, a shift in some suitable feature of the d band structure 
between surface and bulk could provide a suitable generic measure of reactivity. The struc- 
ture of the occupied portion of the d band could be studied experimentally by photoelectron 
spectroscopy, and the shift in the peaks of the spectrum from surface to bulk could be de- 
tected. This measure of the surface d band shift (SDBS) is somewhat ambiguous because of 
broadening and the role of matrix elements. A sharply defined measure is the shift in the 
center of gravity of the d band, B d , or 

SDBS « AB d = -{B s d urf - B b d ulk ) . (1) 

Another sharply defined measure would be the shift in the relevant matrix element of the 
Kohn-Sham Hamiltonian; how to model this is discussed in Sec. III. 

A recent development which is very interesting is that the differences between the core- 
electron binding energies of surface atoms of supported metal monolayers on transition-metal 
substrates and of surface atoms of the clean elemental crystal surface (the latter consisting of 
the element forming the adlayer) are strongly correlated with the corresponding shift in the 
center of gravity of the surface d band [0|§. It has been observed that the core- level shifts 
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tend to correlate with the adlayer's chemical reactivity and that this could be understood 
through the correlation with the surface d band shift 

The SCLS's and the SDBS's arise from the surface potential shift experienced by the core 
states and the d states, respectively. Both represent different spatial samplings of the same 
potential shift; thus the existence of the observed correlation is not surprising. The purpose 
of this paper is to provide a quantitative test of the correlation and a deeper and more 
detailed understanding of its origin. In Sec. II, we report the results of SCLS calculations 
for the 4d transition metals and for Sr and Ag. We also demonstrate an excellent quasilinear 
correlation between the computed SCLS values and those of the SDBS's. Indeed, the SCLS 
values nearly equal the SDBS values in the initial-state approximation. To demonstrate that 
this correlation is robust, holding for other measures of the SPS, we extract values of that 
shift for the xy and x 2 — y 2 subbands at the center of the surface Brillouin zone by fitting 
the computed bands to a simple model in Sec. III. The correlation is found to persist. 

The two most prominent features of the results in Sec. II and III are (i) the quasilinear 
correlation between the SCLS's and the SDBS's and (ii) the sign change which occurs across 
the 4d series. Understanding (i) requires less depth of analysis than understanding (ii). 
Accordingly, we present in Sec. IV plots of the 3d core-orbital density and the 4d valence- 
orbital density overlaid upon the SPS which make clear that the sampling argument referred 
above is valid. On the other hand, to understand (ii) fully requires probing more deeply 
into the origin of the SPS than is done in the usual surface d-band-narrowing argument. 
In Sec. V we provide such an analysis of the origin of the SPS via a decomposition of the 
SPS into its component parts together with a decomposition of the electron-density changes 
Ap(r) responsible for it. In the process, we confirm the essential correctness of the band- 
narrowing argument. The resulting improved understanding of the SPS and Ap(r) is our 
most significant contribution. We discuss our results briefly in Sec. VI. 
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II. CORRELATION BETWEEN SCLS'S AND SDBS 



We have considered the fcc(OOl) surfaces of the 4d transition metals and of Sr and Ag 
and have calculated the SCLS's, A c , for the 3d levels. All calculations have been done for the 
fcc(OOl) structure to facilitate intercomparison and the establishment of trends, irrespective 
of whether fee is the equilibrium structure. 

Using Slater's transition-state concept M] to evaluate total-energy differences, A c can 
be estimated from 

A c ^-[ e r f K = i.5)- e r ik K = i.5)] , (2) 

where e s c mi and ej? ulk denote the Kohn-Sham eigenvalues of a particular core state of a surface 
or bulk atom, and n c is the occupation number of the core-orbital. In the initial-state 
approximation, the SCLS's are given by Eq. (FJ) with n c = 2. 

The electron density, core-electron eigenvalues, and densities of states are calculated 
by means of an efficient surface-Green's-function technique based on the tight-binding, all- 
electron linear muffin-tin orbital theory within the local density approximation (LDA) of 



the exchange-correlation functional in the Ceperley and Alder form [TT[] as parametrized by 
Perdew and Zunger fll2"| . The details of the method have been described elsewhere p3| . |14| . 
The potentials are calculated selfconsistently within the atomic-sphere approximation (ASA) 
in an intermediate region consisting of the surface layer, three substrate layers, and two 
layers of empty spheres simulating the vacuum-sample interface. This intermediate region 
is coupled to the semi-infinite vacuum on one side and to the semi-infinite crystal on the 
other, with frozen potentials. Calculations are performed for sphere radii chosen so as to 
minimize the total energy of the bulk in the fee structure. No relaxation of the topmost 
layers has been assumed; all interatomic distances in the surface layer, as well as between 
surface and substrate layers are assumed to be the same and equal to those in the substrate. 
To describe the transition-state (see Eq. ||), selfconsistent electronic structure calculations 
are performed under the constraint of charge neutrality. The single impurity problem of 
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the localized core-hole with half an electron missing in the 3d shell can be easily treated 
within the surface Green's function formalism. The charge neutrality constraint within this 
scheme results in a nearly fully relaxed final state, i.e., that the valence charge of the excited 
single impurity nearly equals (Z+0.5)e (Z atomic number). In Table | we compare three 
of our calculated shifts with independent first-principles calculations, which include final- 
state effects to some approximation, and with available experimental data. A number of 
other calculations for the most close-packed surfaces fl9j have achieved a similar degree of 
agreement between experiment and theory using similar but not identical computational 
methods. 



We have also calculated the LMTO potential parameters Cd ||20|| , which correlate closely 
with the center of gravity of the d band, B^, for both bulk and surface and have constructed 
the SDBS's from them according to Eq. ([J) (AB d « AC^). The calculated initial-state and 
final-state core-level shifts for the fcc(OOl) surfaces of the 4d transition metals plotted vs. 
these calculated SDBS's are presented in Fig. |I|. Note the near linearity of the correlation 
between the initial SCLS's and the SDBS's. The smaller screening contributions to the 
shifts, which are related to the downward shifts of the core-levels when deoccupied leave the 
overall linear correlation largely unaffected. 



III. CORRELATION OF THE SCLS'S WITH SURFACE POTENTIAL 

PERTURBATIONS 

We have somewhat arbitrarily selected the SDBS as defined by Eq. ([I]) as the generic 
measure of the response of the d electrons to a generic perturbation and have shown that 
it has a beautiful quasilinear correlation with the SCLS's. However, the d band structure 
is complex, with changes in width, shape and centroid position at the surfaces. We now 
examine the relationship between the SCLS's and the values of a quite different measure of 
the SDBS to test whether the above correlation is relevant despite the complexity of the d 
band. 
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In earlier papers we have shown how to extract and display the individual physical effects 
of the presence of the surface on the electronic structure by examining the kip, symmetry- 
, and layer-resolved density of states (DOS) at k|| = for both clean surfaces and 
overlayers |[22|| . For the fcc(OOl) surface, the xy and x 2 — y 2 subbands do not couple with 
any other low-lying bands. The bulk xy and x 2 — y 2 contributions to the DOS are accurately 
fitted by a simple cosine band and are thus represented by the band structure of an infinite 
one-dimensional chain of sites having only nearest-neighbor coupling, t, of a non-degenerate 
level. For the surface layer and substrate layers underneath, the xy and x 2 — y 2 subbands 
DOS at k|| = 0, are accurately represented by that of a semi-infinite chain perturbed by 
potential shifts V\ and V2 at the terminal and penultimate sites. The quantity V\ can 
be interpreted as a model of the surface shift of the diagonal matrix element of the true 
Hamiltonian in a generalized Wannier representation. We take its negative as an alternative 
measure of the SDBS which reflects in quantitative detail the shift in the surface d band 
structure caused by the surface potential shift. The band widths At were fixed by fitting 
to the bulk subbands and the values of V\ and V2 were fixed by fitting the zeroes in the 
DOS in the third layer. The resulting fits for the full energy and layer dependence were in 
general extraordinary as shown in Fig. I3J for the Rh xy subband, for which the fitted and 



the calculated DOS are nearly undistinguishable. In Ref. pl| , the DOS was fitted with V\ 
only. The fits with two parameters are slightly better, but the V\ values are robust, i.e. 
they change little between the one and two parameter fits. The calculated initial-state and 
final-state core-level shifts for the fcc(OOl) surface of the Ad transition metals plotted vs. the 
fitted surface potential shifts, Vi, for the x 2 — y 2 , and xy subbands, are presented in Figs. 
2a and 2b, respectively. Once again we find a quasilinear correlation between the initial- 
state shifts and an independent measure of the SDBS's which is relatively little affected by 
final-state contributions to the core-level shifts, demonstrating that the correlation is indeed 
robust. 
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IV. THE SURFACE POTENTIAL SHIFT AND THE ORIGIN OF THE 

CORRELATION 



It is very interesting that the numerical values of the SCLS's are not merely correlated 
with, but they are nearly equal to those of the potential shifts —Vi(xy), — Vi(x 2 — y 2 ), 
and ACd, and particularly so for the initial-state contributions to the shifts and ACd- To 
understand better this nearly equality we have examined the spherically symmetric part of 
the surface potential shift AV(r), 



together with r 2 \R 3 d(r)\ 2 and r 2 |i? 4 d(r)| 2 , where R^dir) and Ruif) are the radial solutions 
of the Schrodinger equation within the corresponding atomic sphere in the bulk. To first 
order, 



hold, with appropriate normalization. Eq. (4) is expected to be more accurate than Eq. (5) 
since the 3d core-level is much more tightly bound and the LMTO-orbital by itself cannot 
represent the full complexity of the 4d band. 

We show our results in Fig. 4a for Y, representative of the lighter 4d transition metals 
with a positive shift, in Fig. 4b, for Mo, representative of the mid series elements with 
small shifts, and in Fig. 4c for Pd, representative of the heavier Ad transition metals with a 
nearly filled d band and a negative shift, We compare the computed values for the initial- 
state contributions to the shifts and ACd with those estimated via Eqs. (4) and (5) in 
Table JI[ We see that the agreement is better for the initial-state SCLS than for the AC^s, 
as anticipated, but Figs. 4 and Table clearly imply that the sampling argument for 
explaining the correlation is correct. It could not be anticipated in advance of these plots, 
however, that the sampling argument must be correct, considering the complexity of the 



SPS = AV(r) = V surf {r) - V bulk {r) 



(3) 




(4) 



(5) 



S 



spatial dependence of AV(r) and the fact that the 3d and 4c? orbitals are localized in quite 
different regions of space. However, the spread of the Ad orbitals effectively averages over 
the spatial fluctuations of AV{r). 

It should be noted that for radii r larger than half the nearest neighbor separation, 
indicated by a vertical line in both Figs. 4 and 5, the computed quantities become unreli- 
able because of the use of the ASA. The resulting uncertainties do not weaken any of the 
conclusions drawn above from Figs. 4 and below from Figs. 5. 



V. ORIGIN OF THE SURFACE POTENTIAL SHIFT 



Figs. 1, 2a and 2b show the well-known change of sign that was observed to occur in 



the middle of the 5d transition series |S| and predicted for the most close-packed surface of 
the observed crystal structure of the Ad series [[HJ . This fact has usually been qualitatively 
explained by using the decrease in width of the d band at the surface compared to the bulk 
and assuming an approximate conservation of d charge in each layer. Thus, the self-consistent 
potential at the surface changes so as to mantain the d band filling approximately constant, 
and consequently the surface d band shifts relative to the bulk band. This perturbing 
potential acts on the core electrons as well and is repulsive for late and attractive for early 
transition metals. These arguments apply as well to the nominally empty or filled d bands 



of Sr and Ag [p4 |. For Ag, the sp states below and above the Fermi level, Ep, are hybridized 
with the d states below Ep. The d band narrows at the surface, as above. Were the d 
band center to remain unchanged, the sp — d hybridization of the empty states above Ep 
would decrease, and the net d-character f d of the occupied states would increase. Here fd is 
defined as fd = J Ef dErid(E), where Ud{E) is the density of states per atom projected on to 
the LMTO d basis functions, either for a bulk or a surface site. A repulsive potential shift 
moves the d band center upward and conserves the d character. For Sr, the sp states are 
hybridized with states from the empty d band above Ep. When the empty d band narrows, 
fd would decrease if the center of the d band were to remain fixed. An attractive potential 
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shift moves the d band center downward and approximately conserves fa- 
it is important to recognize that in this argument the potential shifts are assumed to 
occur primarily to conserve d character and not strictly to preserve electrical neutrality. 
For the heavier transition metals, it is known |25|] that the outflow of electrons from the 
surface layer to the vacuum, which generates the surface dipole layer, indeed originates 
predominantly from states of s and p character. Despite the qualitative consistency of the 
d band-narrowing argument, we regard it as incomplete and now present a complementary 
but deeper analysis. The new analysis yields a considerably more detailed understanding 
of the surface potential shift and in the process confirms the essential features of the band- 
narrowing argument. 

The Kohn-Sham potential Vtot is conveniently decomposed into three components in 
our computational methods, the Coulomb term V c , which is the Coulomb potential within a 
muffin-tin sphere arising from all charge within that sphere, the Madelung term Vm, which is 
the Coulomb potential arising from all charge external to that sphere and which is constant 
within the sphere, and the exchange-correlation term V xc which is evaluated in the LDA. 
Thus, the change in the total potential is 

AV T ot = AV M + AV C + AV XC (6) 

where the A indicates the difference between surface and bulk quantities. The AVm are 
constants and are listed in Table III for Y, Mo and Pd. The spherical averages of AVm + AV c , 
AV XC , and AVtot are shown in the upper panel of Figs. 5a-c. One sees inmediately that 
AV XC is of significance only in the outer region of the atomic sphere, where it is repulsive. 
This occurs because to lowest order the derivative 



AKc(r) 



Ap(r) (7) 



dp 

diverges as the density goes to zero, and the density is lowest in the outer regions of the 
sphere. However, comparing AV xc (r) in Figs. 5 with the 3c? and 4d orbital densities in Figs. 
4, one sees that AV xc (r) makes little contribution to the SCLS's and makes a comparable 
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negative contribution to the SDBS's in all cases. AVm is positive and AV C is negative in all 
three cases. Moreover, both AV M and AV C vary monotonically in the series Y, Mo, Pd as 
illustrated by the values of AVm(V) and of AV xc (r) at r = and r = r s (r s is the atomic 
sphere radii within the ASA) listed in Table III and by the curves in Figs. 5. It is only the 
sum of the Coulomb and Madelung potential shifts which establishes the trend in the SPS's, 
the SCLS's, and the SDBS's. 

This behavior of the potential shifts is associated with a corresponding selfconsistent 
behavior of the spherically averaged electron density p(r). The shifts Ap(r) are plotted in 
the lower panels of Figs. 5 for Y, Mo and Pd, respectively, together with their decomposition, 
Api(r), into contributions from states of given angular momentun, I = s, p, d. One sees 
that the s and p contributions to Ap(r) are everywhere negative. Table IV lists the electron 
number shifts 

AQ l = 8n r r 2 A Pl {r) dr (8) 
Jo 

and the total shift per sphere 

AQ = J2^Qi (9) 

for the surface layer together with AQ for the first vacuum layer. One sees first that the 
AQ's are essentially equal and of the opposite sign for the vacuum and surface layers. That 
is, the net electron transfer occurs primarily between the surface and the vacuum, with little 
transfer occurring between the surface and the interior. Next, one sees that for Pd 

\AQ S + AQ P \ > \AQ\ (10) 

in the surface layer. Thus, for Pd, s and p electrons both flow out of the surface layer and 
are transferred into the surface d bands, giving rise to the positive AQ d . For Mo, there is 
less than one third as much internal transfer from s and p to d, leading to a small AQd- 
Finally, for Y, d electrons flow out of the surface into the vacuum, giving rise to a negative 
AQ d - 
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This trend in AQd is precisely what is expected from the d band-narrowing argument. 
Because for Pd the band centroid is below E F , the narrowing of the surface d band initiates 
a transfer of s and p electrons into d electrons states at the surface which tends to increase 
the Coulomb repulsion V c and shift the d band centroid upward thus moderating but not 
eliminating the transfer. For Mo, with E F near the centroid, the effect is much smaller. 
For Y, with the centroid above E F , the narrowing causes a decrease of the d electrons, 
reducing V c and shifting the centroid downwards. One cannot at this point distinguish 
between a transfer of d to s and p electrons which then flow outward and a direct outflow of 
d electrons; the effect is the same. However, the large value of r 2 i?| d (r) in the outer region 
of the atomic sphere evident in Figs. 4a for Y suggests at least some direct d outflow. 

Examination of Ap d (r) in the density-shift plots in the lower panels of Figs. 5 reveals 
more interesting results. One sees clearly that both transfer from sp to d states into the 
inner region of the atomic spheres and d flow out of the outer region into the vacuum occur 
in all three cases. In Pd transfer within the interior dominates, in Mo the two effects nearly 
balance, and in Y the outflow dominates. 

VI. SUMMARY AND DISCUSSION 

We have established in Sees. II and III a robust quasilinear correlation between the 
surface core-level shifts and two different measures of the surface potential shift of the entire 
Ad transition series plus Sr and Ag. We then demonstrated in Sec. IV for the representative 
elements Y, Mo and Pd that the initial-state contributions to the SCLS's are accurately 
given by the average of the SPS over the 3d LMTO-orbital and that the SDBS is given 
approximately by the corresponding average over the Ad orbital. The latter, being quite 
broad, averages over the spatial fluctuations in the SPS, and so, in effect, samples a potential 
shift little different from that sampled by the 3d orbital. This provides a detailed explanation 
of the correlation between the SCLS's and the SDBS's and does not assume an approximate 
spatial constancy of AV(r) throughout the region sampled by both orbitals. Finally, we show 
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in Sec. V that the sign change in the shifts across the Ad series is indeed correctly given by 
the standard band-narrowing argument, but that the situation is considerably more complex 
than envisaged in the argument in its original form. In the heavier elements, the s and p 
character of the surface electron density is reduced both by flow into the vacuum and local 
transfer into surface d states. That is, transfer from s and p states provides the increase in 
the d electrons predicted by the band-narrowing argument. In the lighter elements, it is the 
outflow of d electrons into the vacuum which provides the decrease in d electrons predicted 
by the band-narrowing argument. In the middle of the series, there is little net change in 
the number of d electrons because the transfer from s and p electrons into d states at the 
surface is balanced by the outflow into the vacuum It is this systematic variation in the d 
electron density shift and the shift in the total number of d electrons which is responsible 
for the sign change in the SCLS's and the SDBS's through its contribution to the shift in 
the total Hartree potential AVm + AV^. 

In the introduction, we pointed out that the SDBS could be regarded as a response of 
the d band to a surface perturbation, AVtotM, and as such could be used as a measure 
of surface chemical reactivity. However, it is a very crude measure, indicating intensity of 
response but giving no indication of the spatial variation or the geometry of the response. 
Recently, some progress has been made in understanding the geometry of response as well 
as its intensity through the introduction of the concepts of chemical reactivity theory into 
the discussion of the surface chemistry of metals [p^-p8| . The success of the present study of 
trends in the SCLS's and SDBS's suggests that studies of trends in these more sophisticated 
measures of chemical reactivity could be both feasible and fruitful. 

M. V. Ganduglia-Pirovano thanks P. Fulde for his hospitality at the Max Planck Institute 
in Stuttgart, where part of this work has been carried out. One of us (J.K.) acknowledges the 
financial support from the Grant Agency of the Academy of Sciences of the Czech Republic 
(Project No. 110 437) and the Austrian Science Foundation (Project P10231). 
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TABLES 



SCLS (eV) 

Calculated Experiment 
Rh(OOl) -0.83 a ,-0.62 b -0.62 c 

Pd(OOl) -0.34 a ,-0.33 fe -0.44 d 

Ag(OOl) -0.11 a ,-0.10 b -0.0 + 0.1 6 [fcc(lll)] 

TABLE I. Comparison between independent first principles calculations of surface core-level 
shifts with available experimental data. 
a Present work 

fe Methfessel et al. (Ref. H) 
c Borg et al. (Ref. |7|) 
d Nyholm et al. (Ref. @) 
e Andersen et al. (Ref. [(TSj] ) 







order perturbation theory 


Computations 




SDBS 


A initial 


SDBS 


A initial 


Pd 


-0.031 


-0.042 


-0.042 


-0.046 


Mo 


+0.010 


+0.013 


+0.011 


+0.018 


Y 


+0.021 


+0.044 


+0.038 


+0.038 



TABLE II. Comparison of the first order estimates for the 3d coredevel and d band shifts with 
computed initial-state shifts and d band shifts, in Rydbergs. 
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r = 






T = Tr, 

' 1 s 






AVm 


AV C 


AVm + AV C 


AVm 


AV C 


AVm + AV C 


Pd 


+0.18 


-0.16 


+0.02 


+0.18 


-0.14 


+0.04 


Mo 


+0.26 


-0.26 


+0.00 


+0.26 


-0.28 


-0.02 


Y 


+0.16 


-0.18 


-0.02 


+0.16 


-0.21 


-0.05 



TABLE III. Potential shifts in Rydbergs at the center and surface of the muffin-tin sphere. 





s 


Surface 

AQi 
P 


d 


AQ 


Vacuum 
AQ 


Pd 


-0.090 


-0.227 


+0.077 


-0.240 


+0.248 


Mo 


-0.075 


-0.330 


+0.022 


-0.382 


+0.385 


Y 


-0.035 


-0.170 


-0.117 


-0.322 


+0.300 



TABLE IV. Shifts in the number of electrons per atom at the surface and vacuum layers. 
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FIGURES 

FIG. 1. The calculated initial-state and final-state core-level shifts and the the shift of the d 
band center between the bulk and surface d bands. 
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FIG. 2. The calculated initial-, and final-state core-level shifts and the fitted surface potential 
shifts, Vi, for the (a) x 2 — y 2 and (b) xy subbands. 
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FIG. 3. The solid lines are the calculated xy contributions to the k||-, and layer-resolved DOS 
at k|| = for a Rh(OOl) surface. The topmost three sample layers are denoted si, s2, and s3. The 
dotted lines are the local density of states of the terminal, and first two interior neighbors of a 
perturbed semiinfinite chain with nearest neighbor coupling 2t = 0.0980 Ry, and potential shifts 
V\/2t = 0.4984 and V2/2t = 0.0010 on the terminal and penultimate sites. The position of the 
bulk substrate Fermi level is at —0.0765 Ry. 
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FIG. 4. Spherically symmetric part of the surface potential shift and the radial solutions of the 
Schrodinger equation within the corresponding atomic-sphere in the bulk, r 2 \R n i(r)\ 2 , nl = 3d,4d 
for (a) Y, (b) Mo, and (c) Pd. The vertical lines correspond to a radius of half the nearest neighbor 
separation. 
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FIG. 5. Spherically symmetric part of the contributions to the surface potential shift and 
the r 2 Api(r), I = s,p,d of the surface atoms for (a) Y, (b) Mo, and (c) Pd. The vertical lines 
correspond to a radius of half the nearest neighbor separation. 
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